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1 Introduction 



A charged point-particle creates an electromagnetic field strength, the Lienard-Wiechert 
field, that diverges at the position of the particle. This implies that the radiation reaction, 
or self force, experienced by the particle is infinite, and the Lorentz-equation is ill defined. 
The divergent contribution from this equation can be eliminated by a classical infinite 
mass renormalization, to obtain a finite effective equation of motion, the Lorentz-Dirac 
equation [1], that plays a crucial role in classical radiation theory. This equation takes 
the radiation reaction into account and completes the Larmor formula. 

The Lorentz-Dirac effective equation of motion bears several unusual features. It is of 
third order in the time derivative, and it can not be deduced from a standard Lagrangian: 
eventually it must be postulated. The ultimate justification for the equation comes from 
the requirement of local energy-momentum conservation, i.e. there should exist an energy- 
momentum tensor T^" that is a) conserved, and b) admits finite momentum integrals. The 
naive energy-momentum tensor 6'^'^, being quadratic in the field strength F'^'^, diverges 
at the position of the particle like Q^" ~ if R is the distance from the particle, and 

it does therefore not admit finite four-momentum integrals. Said differently, while F^^ is 
a distribution 0^*^ is not a distribution, because the square of a distribution is in general 
not a distribution. Consequently, the four-divergence dfj,Q^'^ does not even make sense. 

The construction of a consistent energy-momentum tensor requires, in some sense, 
to isolate and subtract from Q'^" the singularity present at the position of the particle, 
without modifying the value of Q^'^ in the complement of the particle 's worldline, in com- 
patibility with energy-momentum conservation and Lorentz invariance. More precisely, 
the so obtained renormalized energy-momentum tensor should be conserved, if the parti- 
cle satisfies the Lorentz-Dirac equation of motion. It is clear that such a program can be 
carried out only in the framework of distribution theory. 

Until now the construction outlined here has been realized only for a charged point 
particle in four dimensions, rather recently [2], using a somewhat cumbersome and implicit 
distribution technique. A physically more transparent alternative representation of the 
resulting energy-momentum tensor - again in the framework of distribution theory - has 
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been given in [3], relying on a new Lorentz-invariant regularization scheme, followed by 
a classical renormalization. 

Aim of the present letter is to generalize the new approach of [3] to construct a 
consistent energy-momentum tensor for a system of point like dyons, taking the radiation 
reaction into account, which has not been given before ^. This result completes the 
consistency proof for a classical system of radiating dyons, satisfying generalized duality 
invariant Lorentz-Dirac equations [4], [5], see (2.8). 

The new method illustrated here, due to its manifest Lorentz-invariance at each step, 
appears in particular suitable for extension to a system of strings or branes in arbitrary 
dimensions. Using this method we hope indeed to furnish elsewhere the construction of 
a consistent energy-momentum tensor for a generic radiating extendend object, that is 
still unknown. 

2 Regularized Maxwell— equations 

For simplicity we consider a single dyon with mass m, electric and magnetic charges e 
and g, and worldline the extension of our construction to a system of N dyons 

being straightforward. We denote four-velocity, four-acceleration and derivative of the 

four-acceleration by = dy'^/ds, = du^/ds, = dw^/ds. Introducing a current 
with unit charge as, 

f{x)^ ju^5\x-y)ds, (2.1) 

the electric and magnetic currents are = ej^ and = gj^. The Maxwell equations 
for the dyon become then, 

d^F^" = fe, (2.2) 
d.F'" = fm, (2.3) 

where we indicate the dual of an antisymmetric tensor with F^'^ = ^e^^f"' Fpa-. The 
general solution of (2.2), (2.3) can be written as, 

F^"" ^r'' + e H^"" - g H^"", H'"' = d^A'' - d'^A^, (2.4) 

^The energy-momentum tensor proposed in [4] requires to modify the naive tensor Q^" also in the 
complement of the worldUne. 
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where f^^'^ is a free radiation field, d^f^^ = = duf^" i and Ai^ is a unit Lienard-Wiecliert 
potential in Lorentz-gauge, satisfying UA>^ — j^, d^A^ — 0, 

A'' = — r^- (2.5) 

We write the scalar products as a'^h^ = {ah), and we have defined, 

Ri'{x)=xi' -yf'is). (2.6) 

The proper time appearing in y^^ and in u^^ is the retarded proper time s{x) determined 
from, 

{x-y{s)f^^, x'>y\s). (2.7) 

This means in particular that we have R^Rn — 0, and hence R^ — \R\ = R. 

Since A^^ carries 1/R singularities, the unit field strength H^'^ carries l/R"^ singulari- 
ties near the worldline, and the Lorentz equation for the dyon would be singular. After 
subtraction of the singularity one postulates the following finite duality invariant gener- 
alization of the Lorentz-Dirac equation for a dyon [5], — mu^, w'^ — w^W/^, 
dpi^ e'^ + g'^ ( dw'' 



+ w'u^) + {er + gr)u., (2.^ 



ds 6n \ ds 

that takes the radiation reaction into account. For g — one gets back the Lorentz-Dirac 
equation. 

On the other hand, while A and F have at most 1/R^ three-space integrable singu- 
larities near the worldline of the particle and are distributions, the naive electromagnetic 
energy-momentum tensor ^, 

^ [FFY" - ^ ri'^'iFF), (2.9) 

carries thrcc-spacc non integrable 1/ R'^ singularities, and is not a distribution. 

As first step to isolate the singularities in Q^'^ we introduce a Lorentz-invariant regular- 
ization, parametrized by a positive regulator with the dimension of length £, by replacing 
the retarded proper time s{x) appearing in A^ in (2.5), with a regularized retarded proper 
time Se{x), determined from, 

{x-y{s)f^e\ x'>y\s). (2.10) 

^In the following we use for the contraction of two antisymmetric tensors the notation [ABY" = 
Ai^pBP", and {AB) = Ai^^B^^,. 



We call the resulting regularized potential, 

, a^^^^o, (2.11) 



- 

' 47r(^i?) ^^^^^ 

where from now on with y'^, m'^, R'^ etc. we intend their regularized versions, obtained 
through the replacement s{x) — > Se{x). The regularized field strength becomes then, 

p^u ^ji^^^e - g H^", (2.12) 

where, 

= a'^A^ - a^A^. (2.13) 
We define the regularized energy-momentum tensor as, 

ef = {F.F.r-' - ^ r^iF.F,). (2.14) 

The fields A^, and are now all regular on the particle's worldline, indeed they 
are C°°-functions on R^. But, whereas and F^ for £ ^ tend to A and F in the 
distributional sense, converges to © pointwise away from the worldline, but not in the 
distributional sense, because © is not a distribution ^. 

3 Construction of a finite energy— momentum tensor 

Before taking the distributional limit of ©^ wc must therefore separate and subtract its 
contributions that diverge for e — > in the distributional sense. This means that we have 
to apply ©£ to a generic test function (p, and isolate the terms of the resulting integral that 
diverge for £ — > 0. The rest of this section is mainly devoted to the exphcit identification, 
and subtraction, of these singular terms, the result being formula (3.5). 

To begin with we need an explicit expression for F^, i.e. for H^. Differentiating (2.10) 
to derive d^s^ — Rn/{uR), from (2.11) one obtains, 

= ^^^^^^3 [{r]>'P - u>'uP)RX + {u'^w'' - u''w'')RpR''] . (3.1) 



^Saying that a set of functions, or distributions, /e converges for e — > to / in the distributional sense 
means that it converges if applied to an arbitrary test function, i.e. one has linie^o fs{v) = f{v) for 
every G iS(R^). 
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The most singular terms in H^, see (2.13), go therefore as near the worldUne, and 
since the radiation field f^" is supposed to be regular, the contributions of that do 
not converge for £ ^ in the distributional sense, are only the ones quadratic in H^. 
Inserting (2.12) in (2.14) and keeping only the terms quadratic in H^, one obtains for the 
divergent part of therefore. 



(3.2) 

div 



Notice that the cross terms in e and g canceled. 

We remain then with the evaluation of {H^H^)ciiv This product contains terms that 
behave as 1/R"' near the worldhne, with n = 2,3,4. As £ —> in the distributional 
sense, for dimensional reasons the terms with n = 4 give rise to simple pole [~ 1/^], and 
logarithmic [~ In e] singularities, the ones with n = 3 give rise to logarithmic singularities, 
while the ones with n — 2 are convergent. Actually, it can be seen that the logarithmic 
singularities cancel between the n = 4 and n — 2> terms, and one remains only with the 
pole singularities contained in the 1/ terms. To determine {H^H^)iiiv it is then sufficient 
to evaluate the simple pole term of, see (3.1), 



V ^ ^Jdiv 1fi7r2 ^ ' ^ ^ ' ^ n.R\e 



(3.3) 

l/e 



where we kept only the 1 / R'^-terms. The rest of this section is devoted to the explicit 
evaluation of the r.h.s of this formula. 

Since the divergences as £ — > are intended in the distributional sense, we must 
apply the r.h.s of (3.3) to a test function. Omitting for simplicity of writing the (regular) 
tensorial prefactor (77'^^ — tt'^w^) • • • in (3.3), we have to evaluate the function RpRa/ {uR}^ 
applied to a generic test function (fi{x) e S, i.e. by definition, 

^^^^ " / (^ ^^^^ ^^^^ = / / (^ ' - '^^^^^'^ + ^^-^^ 

where we have inserted a (^-function to take the constraint (2.10) into account, and we 
have performed the shift — > x^ + y^. The x-integration is restricted to x° > 0. The 
kinematical quantities y and u are now evaluated at s, that is a free integration variable. 
The pole singularity of this expression can be extracted by rescaling a;'* —>■ sx^^, and 
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sending in the integral e to zero ^, 



RpRa 

(uR)^ 



i/e — W 

This means that, as expected, the part that diverges as £ — >■ in the distributional sense, 
is entirely supported on the worldline. Indeed, we can write the above result as. 



RpRa 



{uRf 



l/e 



^ j {bUpUa- - r]pa) S^{x - y) ds. 
Inserting this expression in (3.3) we get, 

{d^^A'^d'^A^)^^ " eJi / ^^"^"^ ~ 5\x - y) ds. 

Using this result in (3.2) allows one to determine the divergent part of the energy- 
momentum tensor as, 

©rid.. = / (^'^^^ - \ r)'"") - y) ds. 

This means that we can define a "renormalized" energy-momentum tensor for the 
electromagnetic field as, 

,2 I ^2 



^S' - lim 



■ em 



er - / {u^u'' - \ r") s'{x - y) ds 



(3.5) 



where S' — lim means limit in the distributional sense. What we have shown here is 
that this limit exists and represents a well defined distribution. This means in particular 
that the four-momentum integrals over an arbitrary finite three-volume Jy d^xT^^ exist, 
whether or not the particle at the given instant is inside V. If the acceleration of the par- 
ticle vanishes sufficiently fast for t — > — oo, then also the total four-momentum is finite [2], 
see eq. (4.9) for an expficit expression. Notice also that the expression for T^^, apart from 
being manifestly Lorentz-invariant, coincides with the naive tensor Q^'^ in the comple- 
ment of the worldline. This feature realizes the requirement that the energy-momentum 
tensor should be "changed only at the position of the particle" . The counterterm u'^u''' 
in (3.5) can be interpreted as a kind of mass term, while the term proportional to rj^^ is 
needed to keep T^^ traceless. 



^The integral (3.4) contains actually also logarithmic divergences ~ Ine, but these can be seen to 
cancel against the logarithmic divergences present in the terms, as mentioned above. 
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4 Energy— momentum conservation 

Above we have constructed an energy-momentum tensor that gives rise to finite momen- 
tum integrals. Given the (a priori) arbitrariness of our construction, its physical justifi- 
cation arises from the fulfillment of local energy-momentum conservation. The check of 
this conservation law requires the evaluation of the four-divergence d^T^^. The present 
section is devoted to this evaluation, the result being given in (4.8). 

We begin by stating the form the (regularized) Maxwell equations satisfied by F^'^. 
Prom (2.12) and (2.13) one obtains, 

d.Fr^ej:, d.Pr^gj:, (4.1) 

where the regularized unit current, 

= (4.2) 
which is still conserved, can be calculated applying one more derivative to (3.1), 

The factor of e'^ arises from the fact that R^R^ = e"^ and it implies that for £ ^ 0, in the 
complement of the wordline converges pointwise to zero. More precisely, one has the 
distributional limit S' — lim£_»o if — i^' ^ imphed by (4.2). 

We come now back to the evaluation of d^T^^. Since the convergence in (3.5) is in the 
distributional sense and since the distributional derivative is a continuous operation, we 
can interchange the derivative with the limit, 

d,TZ = s'- Jim (^a,er - / (^^ - i ^'(^ - y) 

Using the regularized Maxwell equations (4.1) one obtains, 

a,ef = -{eF:^ + gF:^)j,, (4.5) 
= -(e^^g')H^/3e,-{er + 9r)3e,. (4.6) 

In the terms containing the external field one can simply take the limit jfj,, whereas 

in the first term one has to evaluate carefully the distributional limit of the product 
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H^^jei^- Due to the factor £^ in front of (4.3), this product converges pointwise to zero 
in the complement of the worldhne; this means that its distributional limit - if it exists 
- is necessarily supported on the worldline. On the other hand, the distributional limit 
5'— lim^^o H^f^j^^ cannot exist. Indeed, since the four-divergence of a distribution is again 
a distribution, we know that the limit in (4.4) exists, and this implies that the divergent 
contributions of H^^j^^ must precisely compensate the 1/e terms in (4.4). Thus, actually 
it is sufficient to determine the finite contributions of H^^j^/^, as £ —> in the distributional 
sense. 

Using the techniques illustrated in the appendix it can indeed be shown that one has, 

1 fdw'' o ,\ 1 / 1 



5%x-y)ds+o{e), (4.7) 



where o{e) stands for terms that go to zero as £ in the distributional sense. Using 
this in (4.6) one sees that the 1/e terms in (4.4) cancel, and one obtains. 



d, 



S'^{x-y)ds. (4.8) 



This implies eventually that, when adding the energy-momentum tensor of the dyon, 
j^nu _ rp^u^ + mj u'^u'^5'^{x — y) ds, the total energy-momentum tensor is conserved, if the 
generalized Lorentz-Dirac equation (2.8) holds. 

If the external field vanishes, f^^ = 0, cq. (4.8) can be integrated over whole three- 
space at fixed time t, to obtain the derivative of the total four momentum of the elec- 
tromagnetic field, If the four-acceleration of the dyon vanishes sufficiently fast for 
t — oo, a further integration gives then the total four-momentum of the electromagnetic 
field as, 

PLit)^Jd'xT^i:,^-^^±f(^w^{s) + £^w\X)u^{^^ (4.9) 
where s is the proper time of the particle at the instant t. 

5 Interpretation and Outlook 

The construction of the energy-momentum tensor performed here supplies further evi- 
dence of the consistency of the classical dynamics of a radiating dyon system. Indeed, for 
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a system of particles our construction generalizes simply by replacing in (3.5) the coun- 
terterm with ^ Y^ri^r + 9r) I ('^r'^^r ~ i V'^"^ ^^(^ ~ Vr) dsr- This subtraction is sufficient 
since the mutual interactions do not give rise to singularities in 0^*^. 

A question that arises naturally is whether the energy-momentum tensor constructed 
in (3.5) is determined uniquely. If we insist on the physical requirement that off the 
worldline this tensor should coincide with the original one (2.9), and if we enforce duality 
invariancc, a priori the electromagnetic energy-momentum tensor is indeed determined 
only modulo a term, supported on the worldline, of the form, 

^TZ = {e' + /) / h'^^5\x - y) ds, 

where h'^" is a symmetric tensor, of dimension one over length, constructed with m^, w^, 
dw^ /ds etc. The question is now if there exists a tensor h^^ for which the modified total 
energy-momentum tensor T>^'^ = + + rn J u^u^5'^{x — y) ds is still conserved. 

Given the above form of ATg'^ the answer to this question is the same as in the case 
of charged particles, and it has been given in [2]: there exists no h^'^ ^ such that 
d/^T^" — 0, whatever modified Lorentz-Dirac equation one imposes on the particle. Given 
the above requirements, and the implicit assumptions that our dyon is point-like and 
spinless, four-momentum conservation fixes therefore (3.5) uniquely. 

Our regularization appears also particularly useful for the derivation of effective equa- 
tions of motion for extended objects, alternative to [6]. In the present case e.g. the self 
force can be obtained evaluating F/*^ in (2.12) at the worldhne x — y(s), using (6.1), and 
then taking £ — > 0. The result is, 

m^iyis)) = — (uV - uV) - — - + o(£). 

' " 87r£ ^ ' Q-Ky ds ds I ^ ' 

Using this in the regularized Lorentz equation for dyons, ^ = {eF^'^ + gF^^^ u^, one 
sees that the divergent part renormalizes the mass, and that the finite part amounts to 
(2.8). 

Our regularization scheme admits a simple interpretation in terms of the retarded 
Green function G{x) = ■^H{x'^)5{x'^) of the Laplacian □ = 9^9'*. It is indeed immediately 
seen that the regularized potential (2.11) is produced by the regularized Green function 
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G^{x) — ^H{x^)S{x'^ — £^), where H is the Heaviside function, according to A^{x) — 
J d^zG^{x — z)j^{z). This regularization extends naturally to arbitrary dimensions since 
in even space-times, D — 2n + A, the Green-function is G = H{x'^)/2TT^~^^{d/dx'^)"-S{x'^), 
while in odd ones, D = 2n + 3, it is G = H{x'^)/27i''+\d/dx^)'^[H{x^)/V^], see [7]. The 
regularized Green function G^ in arbitrary dimensions is then simply obtained operating 
in G the replacement — > — 

Interpreted in this way our method admits then a natural extension to extended objects 
in higher dimensions. For example, for an electric brane in D dimensions, minimally 
coupled to a p-form gauge field B, we can introduce a retarded regularized potential 
- in Lorentz gauge - according to B^(x) — J d^zG^{x — z)j'^'P\z), where the current 
j(^) is the 5-function on the brane, i.e. its Poincare dual. This potential gives rise to 
a field strength (p + l)-form = dB^, that is regular on the brane, and hence to the 
regularized energy-momentum tensor Q^'^ — ^[{F^F^Y'^ — 2(^+1) V^'^iFeFe)]- Following the 
lines of the present paper it should then to be possible to construct a finite and conserved 
energy-momentum tensor for a generic brane - taking its radiation reaction into account - 
providing thus a physical basis for the effective equations of motion postulated previously, 
[6]. Moreover, the knowledge of this tensor should allow also a systematic analysis of the 
energy-momentum loss of an extended object, due to the emitted radiation. We hope to 
report on this construction soon. 

Acknowledgements. This work is supported in part by the European Community's 
Human Potential Programme under contract MRTN-CT-2004-005104, "Constituents, Fun- 
damental Forces and Symmetries of the Universe" . 

6 Appendix: Evaluation of H^^js^ for £ ^ 

From (2.13) and (3.1) one obtains, 

H^^ = ^ [R^u'' + {u^w" - u''wP)RpR'' - (/i ^ i/)] . (6.1) 
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From this and (4.3) one sees that all terms appearing in the product H^^j^^ are schemat- 
ically of the form, 

where G is a tensor constructed with w and 6, all evaluated at Se{x), whose tensorial 
structure we do not indicate explicitly. One has to apply this expression to a test function 
and to consider the limit £ — > 0. Proceeding as in (3.4) one obtains, 

= e'JdsJj^^^^^2 5ix'-e')x'^^--.x^-G{s)^ix + y) (6.3) 
= / d''^ J (^^1^ 2 5{x' -l)x^^...x^- Gis) b(y) + e x'^d^v) + •••], 

where we have rescaled x ^ ex, and expanded ip{ex + y) in powers of e. The values of 
M and N appearing in H^^j^^ are such that M — N — 4,5,6,7. Therefore, as £ — > 
at most the first three terms of the series above give a non vanishing contribution. To 
conclude the evaluation of Ie{^) one must eventually perform the integration over d'^x. 
This integration can be performed by taking multiple derivatives w.r.t u^, considered as 
an independent variable, of the generating function. 



/ 



,4 25{X" - 1} TT 

d X 



3/2 r 



{uxY (M2)n/2 r(^) ' 

and setting eventually — 1. The x-integration gives thus rise to polynomials in u^. 
One sees then that as £ ^ Ie{'^) reduces to a (finite) sum of terms of the kind 
/ L(s) d ■ ■ - d ip{y) ds, where 1 — 0,1,2 and at most two derivatives on ip appear. This 
means that I^, as £ goes to zero in the distributional sense, is supported on the worldline, 
becoming a sum of terms of the type 1/e'' J L(s) d ■ ■ ■dS'^{x — y) ds. More precisely, the 
terms in H^^j^^ with M — N = 4 converge to zero as £ — > 0, the ones with M — N = 5,6 
give rise to finite and simple pole contributions, while the ones with M — N — 7 give rise, 
a priori, also to double pole contributions. However, by direct inspection one sees that 
the double poles cancel. Indeed, from (6.1) and (4.3) one sees that the terms in H^^^j^^ 
with M — N = 7, are given by /J = (4^)2'^(„j^)8 -R/^ iv^'^ ~ u^W-") . According to (6.2) and 
(6.3) we have (M = 8, AT = 1), 
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and the double pole cancels since / -^^2 5{x'^ — 1) = ff ""/i- Only finite and simple 
pole terms survive then in H^^j^^, and a straightforward but a bit lengthy calculation 
gives (4.7). 
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